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Abstract 

The projected shell model analysis is carried out using the triaxial Nils- 
son+BCS basis. It is demonstrated that, for an accurate description of the 
moments of inertia in the transitional region, it is necessary to take the tri- 
axiality into account and perform the three-dimensional angular-momentum 
projection from the triaxial Nilsson+BCS intrinsic wavefunction. 
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The major advancement in the studies of deformed nuclei has been the introduction of 
the Nilsson potential It was shown that the rotational properties of deformed nuclei 
can be described by considering nucleons to move in a deformed potential. The description 
of the deformed nuclei in medium and heavy mass regions is impossible using the standard 
(spherical) shell model approach, despite the recent progress in the computing power. The 
Nilsson model has provided a useful nomenclature for the observed rotational bands. It 
is known that each rotational band is built on an intrinsic Nilsson state. The Nilsson or 
deformed state is defined in the intrinsic frame of reference in which the rotational symmetry 
has been broken and in order to calculate the observable properties, it is required to restore 
the broken symmetry. 

The rotational symmetry can be restored by using the standard angular momentum pro- 
jection operator 0. This method has been used to project out the good angular momentum 
states from the Nilsson+BCS intrinsic state see also the review article || and refer- 

ences cited therein. In this approach, the angular momentum projection is carried out from 
a chosen set of Nilsson+BCS states near the Fermi- surface. The projected states are then 
used to diagonalize a shell model Hamiltonian. This approach referred to as the projected 
shell model (PSM) follows the basic philosophy of the standard shell model approach. The 
only difference is that, in the PSM, the deformed basis is employed rather than the spherical 
basis. This makes the truncation of the many-body basis very efficient, so that the shell 
model calculations even for heavier systems can be easily performed. 

The PSM approach has been used to describe a broad range of nuclear phenomena 



such as backbending J/J, superdeformed |||| and identical bands |10| with considerable 
success. The assumption in the PSM approach has been the axial symmetry for the deformed 
system to keep the computation simple. In fact, this is a reasonable approximation for well 
deformed nuclei. However, for transitional nuclei, this assumption is questionable. The 
inadequacy of the axially symmetric basis has been clearly demonstrated by moments of 
inertia (the backbending plots) of the transitional nuclei in the rare-earth region. It has been 
shown || that, in the low spin region, observed moments of inertia for lighter rare-earth 
nuclei (for instance 156 Er, 158 Er, 158 Yb and 162 Hf) and for the heavier rare-earth nuclei (for 
instance 172 W, 174 W and 176 W) increase quite steeply with increasing rotational frequency 
as compared to the moment of inertia calculated by the axially symmetric PSM approach 
(see Figs. 14-17 in 0). This can be understood by noting that the horizontal (vertical) 
line in the backbending plot represents the rotational (vibrational) limit since the energy 
E(I) as a function of spin I is proportional to I 2 (I). In this sense, the experimental data 
slants towards the vibrational side in comparison with the existing PSM results. On the 
other hand, the spectrum of a triaxial rotor |TT| is known to vary from rotational spectrum 
to a vibrational one as the triaxiality parameter 7 increases from to 30° and, using this 
model, it has been demonstrated (see Fig. 18 in |J) that the backbending plot indeed slants 
towards the vibrational limit when 7 increases. It is therefore expected that, by using the 
triaxial basis in the PSM, the moments of inertia and other properties of the transitional 
nuclei can be described more appropriately. As pointed out in IjJ, the major problem here 
lies rather in the ground state band. This part of the spectrum is quite insensitive to 
the configuration mixing since the energy and spin values are still low (J < 10), so that an 
improvement of the ground state by allowing some triaxiality is in order. The purpose of the 
present work is to develop a triaxial projected shell model (referred to as TPSM hereafter) 
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approach for the description of transitional nuclei. This requires a three-dimensional angular 
momentum projection and has not been attempted so far except for a short investigation in 
early eighties [0. We have carried out the three-dimensional projection and would like to 
report our preliminary results. 

The shell model Hamiltonian employed is identical to the one used in the axially sym- 
metric PSM approach ||. It consists of Q ■ Q + monopole pairing + quadrupole pairing 
forces 

H = H -^Y: QlQ, - G M F«P -G Q Y: (1) 

Here, H Q is the the spherical harmonic-oscillator single-particle Hamiltonian with a proper 
I ■ s-force while the operators Q and P are defined as 

$m = £<W£<*, ^ = ^£44 Pt = \Y,Q^cic^ (2) 

a/3 Z a Z a/3 

where the quadrupole matrix-elements are given by 

Q»« a i = 5 NN ,(N 3 m\Q^\N']'m'). (3) 

In Eq. (2), a = {Njm} while a represents the time-reversed state of a. The Hartree- 
Fock-Bogoliubov (HFB) approximation of the shell model Hamiltonian Eq. (1) leads to 
the quadrupole mean-field which is similar to the Nilsson potential. Therefore, instead of 
performing the HFB variational analysis of the Hamiltonian in Eq. (1), the Nilsson potential 
can be directly used to obtain the deformed basis. In the present work, we use the triaxial 
Nilsson potential specified by the deformation parameters e and e' 

H N = H - hu(eQ + A 2 J-^ 2 ), (4) 

to generate the deformed single-particle wavef unctions. It can be easily seen that the rotation 
operator e~ l 2 Jz transforms the Nilsson Hamiltonian Hn into the opposite triaxiality (e' — > 
— e') leaving the eigenvalues unchanged. Later, it will be shown that the projected energy 
is independent of the sign of e' and it is sufficient to consider only the non- negative e'. The 
volume conservation also restricts the range of e and e' values to 

-3<e<! |e'| <V3(1 + |). (5) 

The triaxial Nilsson potential has been solved for the rare-earth region with three major 
shells iV = 4, 5, 6 (3, 4, 5) for neutrons (protons). 

In the next step, the monopole pairing Hamiltonian is treated based on the triaxial 
Nilsson basis. We use the standard strengths for the pairing interaction of the form 

G M = (G lT G 2 ^-^)i, (6) 

where — (+) is for neutrons (protons) while G\ and G2 are chosen respectively as 21.24 and 
13.86 MeV in the rare-earth region. The pairing correlations are treated by using the usual 
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BCS approximation to establish the Nilsson+BCS basis. The three-dimensional angular 
momentum projection is then carried out on the quasiparticle states obtained in this way. 
The three-dimensional angular momentum projection operator is given by 

pLk = J dnD{ IK (tt)R(n), (7) 

R(Q) = e - ia ^e~ ll3j ye-^^ being the rotation operator and D T MK (n) =< vIM\R{Q)\vIK >* 
its irreducible representation where {\vIM >} is a complete set of states for the specified 
angular momentum quantum number IM. Since the spectral representation of the projection 
operator Eq. (7) is represented by 

P I MK = Y.\v IM ><v IK l (8) 

v 

it is easy to see that |$' >= e _i ^ Jz |$ >, i.e. the state of the opposite triaxiality to a state 
|$ >, is projected to give 

PUl®' >= HtK<^ 3 '\* >= {rY 1 ** Hik\* > ■ (9) 

This state differs only by a phase factor from P[ IK \& > and thus represents the same physical 
state. It therefore proves that the result of the angular momentum projection should be 
independent of the sign of e'. We have used this property to check the programming since 
it is a non-trivial relation. Note that this justifies the above-mentioned restriction e' > 0. 
Details of the projection technique and algorithm are discussed in an Appendix of ||. 

In the present work, we have diagonalized the Hamiltonian Eq. (1) within the space 
spanned by {P{ IK \^ >} where |$ > is the (triaxial) quasiparticle vacuum state. The TPSM 
eigenvalue equation with the eigenvalue E 1 for a given spin / thus becomes 

E { H KK> ~ E'N 1 ^) F* K , = (10) 

K' 

where the matrix elements are defined by 

Hjck' =< ®\HPkK>\$ >, N J KK , =< ^P'kk,^ > . (11) 

This TPSM equation has been solved for a range of nuclei in the rare-earth region and the 
results of a selected few are presented in Figs. 1 and 2. 

The deformation parameters e used in Fig. 1 are exactly the same as those used in the 
earlier calculations with the axially symmetric basis ||, i.e. e = 0.20, 0.20 and 0.225 for 
156 Er, 158 Yb and 176 W, respectively. The results with e' = 0.0 in Figs. 1 and 2 represent 
these axially symmetric calculations. The experimental moments of inertia (represented by 
circles) increase very steeply. The calculations with e' = 0.0 on the other hand depict a 
very slow increase and is typical of an axially deformed rotational band. The moments of 
inertia in Fig. 1 become steeper with increasing value of e' and the value close to e' = 0.15 
reproduces the experimental data. Roughly speaking, this e' value corresponds to 7 = 35°. It 
should be noted that the experimental moment of inertia shown in Fig. 1 slightly increases, 
in particular for 176 W, at the higher end, whereas the theoretical moment of inertia shows 
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a drop. This increase in the observed moment of inertia can be explained by noting that, 
at around spin I=12 + , a 2-quasiparticle band (i.e. the s-band) will cross with the ground 
band and the energy of the higher spin states will be depressed, so that the moment of 
inertia will effectively increase. In the present calculations, the projection has been carried 
out only from the ground (i.e. the O-quasiparticle) band and this effect is not taken into 
account. The projection from 2- and higher-quasiparticle states requires further work and 
will be reported elsewhere. 

Fig. 2 shows the moments of inertia for some Os-isotopes. It is known that these isotopes 
are 7 soft with very low-lying 7 bands. It is clear from Fig. 2 that, for 184 0s, the moment 
of inertia is well reproduced with e' = 0.15. For 186 Os and 188 Os, the experimental moment 
of inertia can be explained with e' between 0.10 and 0.15. 

In summary, it has been clearly shown in the present work that three-dimensional an- 
gular momentum projection from triaxial Nilsson+BCS deformed intrinsic wavefunction is 
essential for an accurate description of the transitional nuclei. The moments of inertia of 
these transitional nuclei depict a steep increase as a function of rotational frequency in the 
low spin region and this can only be explained with triaxial deformation of 7 ~ 30° since 
the inclusion of 2- and higher-quasiparticle bands will affect little in this spin region. 

We would like to mention that the present work has been exploratory. For a detailed 
study, the energy-surface needs to be analyzed as a function of e and e' to look for the optimal 
deformation. In the present work, the deformation parameter e has been taken from the 
earlier studies in which the axial symmetry was assumed. In a more consistent treatment, 



both e and e' have to be varied in order to search the energy minimum for spin / = |L2 
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FIGURES 



FIG. 1. The experimental and the calculated moments of inertia (6) are plotted as a function 
of the rotational frequency (u>) for 156 Er, 158 Yb and 176 W. It is very evident that the Expt. moment 
of inertia is reproduced with the triaxiality e' ~ 0.15 for all the three cases. 

FIG. 2. The experimental and the calculated moments of inertia (0) are plotted as a function 
of the rotational frequency (uj) for 184 0s, 186 0s and 188 0s. For 184 0s, the Expt. moment of inertia 
is very well reproduced with e' ~ 0.15. For 186 0s and 188 0s, it can be described with e' between 
0.10 and 0.15. 



7 



70.0 



60.0 - 



© 
C\J 



50.0 



40.0 



o 



1840s 




o Expt. 

• 8=0.0 

-o 8=0.05 

+ 8=0.10 

^8=0.15 



8=0.20 



1860s 




"♦8=0.10 

-±8=0.15 



^ 8=0.20 



□- 
♦- 



0.025 



0.075 



0.025 



0.125 0.0250 



1880s 




60.0 



-o Expt. 

• 8=0.0 

-o 8=0.05 

-♦8=0.10 

^ e'=0.15 
^ 8=0.20 



0.125 



50.0 



40.0 



30.0 



CO 



